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ABSTRACT 

We searched for coherent X-ray pulsations from Aql X-1 in a series of RXTE 
observations taken shortly after a recent outburst. During the course of these 
observations, Aql X-1 passes through an apparent "propeller" phase as its lumi- 
nosity fades to its quiescent value. No pulsations were detected, and we place 
upper limits (ranging from 0.52% to 9.0%) on the fractional RMS amplitude of 
any periodic signal contained in the various data sets searched. This result has 
implications for the geometry of the system, if the quiescent luminosity is due 
to continued low-level accretion. Alternatively, our result supports the idea that 
the quiescent luminosity may be due to thermal emission. 

Subject headings: stars: individual (Aql X-1) — stars: neutron 

1. INTRODUCTION 

The origin of the quiescent emission observed from low-magnetic field transient neutron 
stars (NSs) has posed something of a puzzle. When first discovered observationally (in 
the transient Cen X-4; Van Paradijs et al. 1987), the X-ray emission, fit with a blackbody 
spectrum, was too faint to have originated over the entire surface of a 10 km NS. It was 
suggested that the emission was due to accretion over a fraction (~ few%) of the NS 
surface. As additional quiescent NS X-ray spectra were observed and fit with a blackbody 
spectrum (Verbunt et al. 1994; Asai et al. 1996a,b), the small implied emission areas seemed 
to confirm this view. 

The cause of accretion anisotropy over the surface of a low B-field NS has been 
purported to be the "propeller effect" (lUarionov & Sunyaev 1975; Stella et al. 1986), 
in which at low accretion rates, the magnetic field increases in importance relative to 
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the gravitational field in determining the accretion geometry, perhaps expelling much of 
the accretion from the system when the effect of the magnetic field is comparable to or 
stronger than that of the gravitational field at the radius where the Keplerian frequency 
is equal to the spin frequency of the NS. When this occurs, some fraction of the accretion 
may follow the magnetic field hues to the magnetic poles which, if these are offset from 
the rotational poles, could conceivably produce X-ray pulsations. The efficiency of such a 
propeller has been estimated recently (Menou et al. 1999). However, some observations of 
decreasing torque with increasing luminosity may indicate deviations from this standard 
magnetic accretion scenario in low mass X-ray binaries (Chakrabarty et al. 1997). If, 
after an outburst, accretion continues into quiescence then pulsations might be expected, 
particularly if the NS can be shown to have gone through a "propeller phase" , in which the 
accretion geometry in the vicinity of the NS is dominated by the magnetic field, rather than 
by the gravitational field. No pulsations from these objects in quiescence have yet been 
reported, although they have been searched for (Verbunt et al. 1994; Asai et al. 1996b). 

Thus, until recently, the view of the quiescent luminosity of low-B transient NSs has 
been as due to continued accretion over a fraction of the surface of the NS, perhaps caused 
by modified accretion geometry due to the effects of the magnetic field, which may therefore 
give rise to pulsations in the X-ray intensity. 

A different view of the origin of the quiescent luminosity has recently been described 
(Brown et al. 1998), in which nuclear reactions at the base of the NS crust keep the NS 
core heated to temperatures (~ 10*^ K) sufficient to explain a large fraction, if not all, of 
the quiescent luminosity of these objects (10^^~^^ erg s~^); in this manner, the quiescent 
luminosity can be produced in the absence of active accretion. The discrepancy in the 
emission area can be explained as due to an incorrect blackbody assumption - specifically, 
if accretion is shut-off, metals stratify in the NS atmosphere (Alcock & lUarionov 1980; 
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Romani 1987), producing a pure H photosphere, in which the free- free opacity, which 
decreases with increasing frequency, dominates, permitting higher E photons to originate at 
greater depth, thus at higher temperatures (Rajagopal & Romani 1996; Zavhn et al. 1996). 
This produces an emergent spectrum which is spectrally harder than a simple blackbody, 
and - when described as a simple blackbody - requires a higher T, and subsequently lower 
emission area, to parameterize it. Spectral fits with H atmosphere spectra show that the 
observed X-ray spectra from the three such objects for which data exists (Cen X-4, Aql X-1, 
and 4U 1608-522) imply emission areas consistent with a 10 km NS (Rutledge et al. 1999). 
Thus, both their quiescent luminosities and their emergent X-ray spectra can be explained 
in this manner, without invoking active accretion during quiescence. An attractive feature 
of this explanation is that it accounts for the similar luminosities of these three objects, 
without requiring their quiescent mass accretion rates (a product of the propeller efficiency 
and mass accretion rate through the accretion disk) to be serendipitously similar. 

It remains an open question if accretion (at a level of M<10~^'^Mq yr~^) occurs at all, 
contributing some fraction of the quiescent X-ray luminosity. There are few observational 
means by which this may be investigated. One such possibility is to search for metal lines in 
the quiescent X-ray spectrum, which, if present, would imply that metals are being fed into 
the atmosphere at a rate greater than gravity can stratify it (Bildsten et al. 1992, see Zavlin 
et al. 1996 for X-ray spectra of metallic atmospheres). Such spectra can be investigated 
by the latest generation of X-ray spectroscopy missions (AXAF and XMM) . A second 
method is to search for intensity variability during quiescence. In particular, as described 
above, if pulsational variability at a particular frequency were to be found in quiescence, 
this would support the scenario of active accretion along the magnetic field lines to the NS 
surface. Alternatively, if the NS magnetic field were strong enough to significantly affect the 
photospheric opacity (> 10^° G; cf. Zavhn et al. 1995), this would produce apparent "hot 
spots" near the magnetic poles of the NS which, if offset from the rotation axis, could give 
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rise to pulsations; however, it is thought that type-I X-ray bursting sources have magnetic 
fields below this value. Clearly discovery of pulsations or constraints on the pulsed fraction 
during X-ray quiescence (following the onset of a "propeller") would provide useful limits 
for models of quiescent NS emission. 

The system Aql X-1 lends itself to such a search for pulsations. During the X-ray 
decline of a recent outburst, the intensity at first decreased by one order of magnitude 
over 17 days, then abruptly (beginning at ~ 10^^ erg s~^) by 3 orders of magnitude 
over 3 days, which was coupled with the sudden spectral hardening of the X-ray intensity 
(Campana et al. 1998; Zhang et al. 1998a). The X-ray intensity stopped dropping at a 
luminosity ~ 10^^ erg s~^, comparable to the quiescent luminosity observed previously 
(Verbunt et al. 1994). This was interpreted as evidence for a propcllcr-phase, though 
this remains to be confirmed through repeated observations of Aql X-1 in outburst. If 
indeed this behavior does indicate a propeller phase in Aql X-1, it offers the opportunity to 
search for pulsations during the period when the NS magnetic field substantially alters the 
accretion geometry. 

Here we have searched for, and not found, pulsations in Aql X-1 during and immediately 
following its supposed propeller phase. In Sec. 2, we describe our search method and the 
results of our analysis. In Sec. 3 we discuss these results and our conclusions. 

2. ANALYSES 

2.1. Search Methodology 

We used high time resolution data taken with RXTE /PC A (Swank et al. 1996). The 
PCA detectors have a total geometric area of 6500 cm^, and a nominal energy range of 2-60 
keV. We selected for analysis four observations, based on the results of Zhang et al. (1998b). 
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These were listed by Zhang et al. as observations 8, 10, 11, plus a later observation, 12, 
which was not analyzed by Zhang due to the faintness of the source - which made their 
spectral analysis infeasible, but can still be useful in our search for pulsations. Based on 
Zhang's results, the beginning of the propeller-phase appears to be between observations 9 
and 10, and so our timing analyses of observations 10-12 are expected to reveal pulsations 
created by the magnetic-field modified accretion geometry (if any), while observation 8 is 
included simply for comparison purposes while the source was moderately brighter. We 
list details of the observations in Table 1. Data were obtained from the XTE archive at 
HEASARC. 

For our data analysis, we used the standard FTOOLS/XTE v4.0. Using standard-2 
formatted data, we extracted the X-ray spectrum from each observation and subtracted 
background counts (estimated using pcabackest tool), and examined the counts spectra. 
Based on the spectrum of the faintest observation (number 12), which becomes background 
dominated above ~ 12 keV, we selected for our analysis PHA bins 0-20, corresponding in 
energy roughly to 2-12 keV. We then extracted individual photon events from data in El 25 
mode (which has time resolution of ~ 125/xs), and corrected the photon arrival times (time 
of arrival; TOAs) to the solar system barycenter. 

Zhang et al. (1998b) reported strong evidence for a 549 Hz pulsation frequency, 
observed during ~ 10 seconds following a type-I X-ray burst from Aql X-1. (It is possible 
that this may be the second harmonic of the NS spin frequency.) Similar pulsations have 
been observed from type-I bursts of 5 other sources (1728-34, Strohmayer et al. 1996; 
1636-53, Zhang et al. 1996; 1731-260, Smith et al. 1997; a galactic center source, possibly 
MXB 1743-29, Strohmayer et al. 1997; and 1702-43, Markwardt et al. 1999). Rather 
than restrict our efforts to a small range of frequencies around the previously measured 
value, we chose to do a more general search for pulsations over the broad frequency range 
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0.5 — 1024 Hz. This range includes all previously reported rotation frequencies for type-I 
X-ray bursters. We were able to do this more general search without a significant loss of 
sensitivity for reasons which we explain below, in section 2.4. 

The Aql X-1 system consists of a NS and a companion star, both orbiting the system's 
center of mass with a period of about 19 hours (see section 2.2 below). As a result of this 
orbital motion, the NS is accelerating along our hue of sight. Any pulsations will therefore 
be observed with a doppler shifted frequency which is not constant in time. In a standard 
fast Fourier transform (FFT) analysis, the spectral power resulting from such a signal may 
be spread out over many frequency bins, drastically decreasing the probability of detection. 
For example, consider a 550 Hz pulsar in a 19 hour circular orbit with a maximum 
projected velocity of 100 km s~^. If we calculate a coherent power spectrum using T = 2048 
seconds of data, the power from this pulsar will be spread over as many as 70 independent 
frequency bins (for the worst case orbital phase, the observed pulsar frequency covers a 
range A/ = 0.034 Hz, giving TA/ = 70 bins). To counteract this effect, we attempt to 
remove the acceleration of the pulsar signal in the time domain before performing the FFT. 
Because we do not have complete knowledge of Aql X-l's orbital parameters, we must cover 
the possible orbital parameter space with a number of acceleration trials and repeat the 
Fourier analysis for each. 

For our purposes, the orbit (assumed circular) can be characterized by three 
independent parameters — the orbital period Porb, the projected circular velocity v, and 
the orbital phase Xq of the NS at the start of the observation in question. The projected 
circular velocity is the magnitude of the NS's circular velocity projected along our line of 
sight, V — I'Nssini where i is the angle between the line of sight and the normal to the 
plane of the orbit. Note that we could equally well have chosen the projected orbital radius 
(asini = vPorb/STr) in lieu of this parameter. The orbital phase Xq is measured in cycles. 
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and is therefore in the range [0, 1], with xq — corresponding to longitude = (where the 
observed pulsar frequency would be measured at its minimum). The signal from a pulsar in 
such an orbit would be observed with a doppler shifted frequency 



where /o is the rest frequency of the pulsar, c is the speed of light, and t is the elapsed time 
since phase xq- 

Pulsar orbital acceleration searches are typically carried out (see, for example, Anderson 



is applicable when the data cover only a small part of the pulsar's orbit, or the pulsar is 
very strong. In the first case, the short span of orbit is well approximated by a constant 
acceleration, while in the second case, some spreading of the feature in the power spectrum 
can be tolerated without the signal disappearing into the noise. For the Aql X-1 searches 
reported here, we assumed that neither of these conditions was satisfied. Our acceleration 
searches therefore fully correct for an assumed circular orbit (Porb) ^, 2;o). As compared with 
the constant acceleration method, this circular acceleration method may require more trial 
accelerations to cover a given orbital phase space, but the detection significance is greatly 
increased. This is because the circular method recovers significantly more signal power in 
a single power spectrum bin, more than making up for the increased number of trials (the 
detection significance is exponential in recovered power, but only linear in the number of 
trials) . 

After preparing an energy selected, barycentered photon TOA list, we proceeded with 
the search method as follows. First, we assumed particular values for the orbital parameters 
(i.e., a trial acceleration) from within the search phase space (of orbital period, velocity, 
and initial phase; see Table 2). We corrected the TO As for the assumed acceleration by 
introducing a corrected time i, which is a function of the original time t, such that the 





et al. 1990) by correcting data with an assumed constant acceleration f{t) = fo + ft- This 
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frequency as a function of t is constant. Equivalently, we require the integrated phase to be 
hnear in t: 



Integrating equation (1) for a given acceleration trial {Porh,v,Xo), we see that the i^^ TOA 
ti is corrected to: 



where {3 = v/c and u) = 27r/Porb- 

These corrected TOAs were used to construct a time series which was then FFTed 
and used to produce an estimate of the power density spectrum (PDS). The PDS was 
searched for candidates (frequency bins containing statistically significant excess power). 
The process was repeated for each acceleration trial. The spacing of acceleration trials in 
(orbital period, velocity, and phasc)-parameter space was chosen to allow sensitivity to as 
weak a signal as possible while keeping the computational requirements reasonable. 



Similar, but inconsistent, values for the orbital period of Aql X-1 have been published 
by two groups. Based on observations of Aql X-1 during outburst. Chevalier & Ilovaisky 
(1991) measured the orbital period to be 18.97±0.02 hr. Shahbaz et al. (1998) determined 
the orbital period in quiescence to be 19.30±0.05 hr. In a more recent analysis. Chevalier 
& Ilovaisky (1998) report the period to be 18.9479±0.0002 hr, again measured during 
outburst, but they also report a quiescent period within 0.02% of this value. To be safe, 
we chose to cover a range of orbital periods that encompasses a 3a range in all of these 
measurements (as well as the periods in between): 




(2) 




(3) 



2.2. Determination of Searched Parameter Space 



18.91 < Porb < 19-45 hr. 



(4) 
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Fortunately, we were able to cover this entire range with one trial value of Porb (see 
section 2.3 below). 

The projected circular velocity v has not been measured. We can determine an 
upper limit on v for our search by assuming a value for the NS mass, m-^s = 1-4 
( Mq = one solar mass) , and choosing a maximum companion mass to which we will be 
sensitive. We can then use our knowledge of the orbital period to calculate the orbital 
velocity of the NS. Recently, the true optical counterpart of Aql X-1 has been identified as 
a late K type star (Chevalier et al. 1999). Again, to be safe we decided to cover a range of 
velocities corresponding to a companion star mass as high as 1.0 Mq (for all possible orbital 
inclination angles) . This results in a search range of 

< t; < 130 km s-^ (5) 

Note that our search was also sensitive to larger companion masses in a restricted range of 
inclination angles; e.g. our search covered companion masses up to 2.0 Mq for 0° < i < 40°. 

For the starting orbital phase for each observation, rather than rely on a particular 
ephemeris, we simply chose to search the full range 

< xo < 1 cycle. (6) 

Our search phase space is summarized in Table 2. 

2.3. Seeirches Performed 

The RXTE observations included in this analysis were roughly 9-13 kiloseconds in 
duration. The TOA data in each observation are broken up into two or three cohesive 
sections, separated by gaps due to earth occultations. Each continuous section is at 
least 2048 seconds long, and we chose to study these continuous blocks individually. 
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Observation 10, with one occultation drop-out, is therefore divided into two sections, which 
we call observations 10a and 10b. Observations 11 and 12 are similarly divided (into 11a, 
lib, 11c, 12a, and 12b). We analyzed only the first section of observation 8, for a total of 
eight separate data sets. The photons from a given data set were binned into a 2^^ point 
time series, with a time resolution of At — 244 //s. Over the course of the observations, the 
source count rate decreases from about 670 counts s^^ at the beginning of observation 8 
to less than 10 counts s~^ by the end of observation 12, with a background rate of about 
31-37 counts s~^ (see Table 3). 

Our search method was described above in section 2.1. We now describe in detail the 
method used to determine the grid of trial values for the orbital parameters. Since we must 
discretize a continuous phase space, there will likely be some offset between a signal's actual 
parameters and our nearest trial values. The effect of this offset will be some spreading 
out of the signal power in the power spectrum, since the time dependence of the signal's 
frequency will not have been completely removed. In choosing the orbital trials, the general 
idea was to space them just finely enough to keep the remaining frequency drift within 
tolerable limits (to be quantified below). 

We began analytically, and then made some empirical adjustments. Considering 
each of the three orbital parameters {PoThjVjXo) separately, we took partial derivatives of 
equation (1), to determine the frequency drift that results from offsets in the individual 
parameters. For example, an error in the velocity {dv = cd(3) causes a drift of 

df = /o cos (out + 27rxo) d(3. (7) 

We are interested only in the extrema of the frequency range covered by the signal over 
the course of the observation time T, and for simplicity, we chose to make our velocity 
trial values independent of the other parameters. We therefore replace the cosine factor in 
equation (7) by its maximum change over the course of an observation (T = 2048 s), and 
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we replace /o by our maximum search frequency (/o — > 1024 Hz). The maximum frequency 
drift caused by a velocity offset is therefore 

9/max = 1024 0.19 dl3 (Hz). (8) 

A "tolerable" drift must be no greater than the independent frequency resolution of the 
power spectrum, F. Since we are considering each orbital parameter separately, for these 
initial calculations we restrict the error-induced drift from each to F/2. This is somewhat 
arbitrary, and the actual spacings used were decided on with the help of simulations. Note 
that the allowable spacing of the trials can be twice the maximum allowable offset. We now 
have for the velocity trial spacing 

cF / F \ 

Avtriai < = 0.77 . km s'^ (9) 

trial ^ ^Q24 0.19 \5x 10-4 Hz/ ^ ^ 

For a coherent FFT of 2048 seconds of data, the independent Fourier resolution is 
F — 1/T — 4.88 X 10"^ Hz. The search strategy that we settled on actually used a larger 
frequency resolution (see below); for now we will leave the expressions for the trial spacings 
in terms of F explicitly. 

Again, for simplicity, we chose to keep the spacing of the Porb and xo trials independent 
of Porb and Xo, but we did allow for velocity dependence. For the orbital period, we 
ultimately find 



(sTIFms) (looi^)" 



(10) 



Even for the smallest Fourier resolution (P = 4.88 x 10-^ Hz) and the largest velocity in 
our search range {v — 130 km s"^), the orbital period trial spacing turns out to be larger 
than our target search range of (19.45 hr — 18.91 hr) 3600 s hr""^ = 1944 s. In other words, 
a single trial value was sufficient to cover our entire search range in Porb- Thus, we did 
not actually search over trial values of the orbital period. For the initial orbital phase trial 
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spacing, we find 



Axtriai < 1-2 X 10-^ 



5 X 10-4 Hz 



F 




100 km s 



V 



cycles. 



(11) 



If we were to use these spacings with coherent 2048 second FFTs, each of the eight 
data sets would require over 95 000 trial accelerations. This would require > 100 days of 
CPU time on a Sun Ultra 10 workstation. To reduce the computational requirements, we 
chose to utihze incoherently stacked power spectra. The original data of duration T are 
divided into S sections of duration T/S. Each section is FFTed individually and the S 
individual power spectra are (incoherently) added together. The independent Fourier step 
size has been increased by a factor oi S {F — 1/{T/S) — S/T) and the trial parameter 
spacings increase by the same factor. The total number of trials is therefore reduced by a 
factor of 5^ (since we are searching over two parameters) . The reduction in the number of 
trials comes at the cost of reduced sensitivity, so the number of stacks S should be kept 
as small as possible, just barely bringing the number of trials to a feasible level. For our 
analyses, we chose to use 5 = 4. 

With the above analytical calculations as guides, we used simulations to determine 
the trial spacings actually used in the search. Our simulations upheld the decision to use 
a single Porb trial. Since there are indications that the orbital period is more likely to lie 
towards the lower end of our search range (18.91 — 19.45 hr), we chose to use Porb — 19.00 hr 
in our trials. We decided to use a velocity trial spacing of 



This is about 1.5 x finer than the spacing calculated above (recall that we are using S — 4: 
stacks). For the orbital phase, we decided on a trial spacing of 



= 2.0 km s-^ 



(12) 




(13) 
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which is about 2.3 x coarser than the spacing indicated in the analytical calculation. For 
small velocities, we used a minimum of 16 Xq trials (except for the single v — trial). 

Using these spacings with 2x oversampled power spectra (exactly like our actual 
searches), we determined that at least 98% of the time, we were able to recover at least 
77% of a simulated signal's power in a single spectral bin, even in the least favorable 
regions of our search phase space. We simulated signals whose orbital parameters were 
offset from our search trials and whose pulsation frequencies were offset from our discrete 
Fourier frequencies. The reduction in recovered signal power is due to a combination of 
these factors. (Note that it is mere coincidence that a simple FFT recovers, on average, 
77% of a signal's power based solely on offset from the discrete Fourier frequencies. In a 
2x oversampled spectrum, the minimum power recovered, based solely on frequency offset, 
is 81%.) Our peak detected power was not always in the frequency bin closest to the rest 
frequency of the simulated pulsar. This has no effect on the detection of pulsations; the 
true frequency and orbit of the pulsar can be refined after the initial detection. 

In total, we covered the orbital phase space with 4069 acceleration trials. The same 
accelerations were used for each of the eight data sets. 

2.4. Estimation of Detection Sensitivities 

To characterize the sensitivity of our search, we wish to place quantitative limits on 
the minimum signal strength required for a significant detection. Since the noise statistics 
of our power spectra are well understood, it is a simple matter to determine the detection 
threshold, the minimum spectral power Pjet required for a detection. This information can 
be used to determine the minimum required signal strength. For a review of detection 
thresholds, detection sensitivities, and upper limits, see Vaughan et al. (1994). 
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Our stacked power spectra were constructed by summing four individual spectra, each 
normalized to a mean power of one (note that this convention differs from the commonly 
used Leahy normalization (Leahy et al. 1983), for which the mean noise power is two). 
In the absence of a signal, the power P in a given spectral bin follows a distribution. 
Specifically, 2P is distributed with IS — % degrees of freedom. The probability that the 
power P in a single PDS bin will exceed a given value Pq is therefore 

p(P > Po) = e-^° (^1 + Po + \PI + \pi^ . (14) 

Considering each frequency bin of each PDS to be one search "trial," our entire 
search consisted of A^triais — 3.4 x 10^° such trials. This is the product of the number of 
observations searched (8), the number of acceleration trials for each observation (4069), 
and the number of frequencies searched in each power spectrum (2^°). The latter number 
is equal to the frequency range searched (1024 Hz) divided by the independent frequency 
resolution {F — (5/2048) Hz) times the oversamphng factor (2). Of course, neighboring 
bins in an oversampled spectrum are not independent, and neighboring acceleration trials 
may not produce truly independent spectra. Thus iVtriais is an upper limit to the number 
of statistically independent trials in the search. Since we are overestimating the number of 
independent trials, our detection threshold will be conservative. 

The statistical significance of a measured power Pq is equal to the probability that 
the power was produced by a random noise fiuctuation. For small p (large powers), the 
significance is 

^ = iVtnalsP(P>Po). (15) 

To achieve our target significance of 10~^, we require a PDS power of at least 

Pdet = 43. (16) 

For a given data set containing total counts, Ng of which were emitted by the 
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source, we can relate the source strength to the expected spectral power by 

1 Ar2 

(Buccheri et al. 1987; Van der Klis 1989) where T is the pulsed fraction, i.e. the fraction 
of source counts that actually contribute to the pulsation. The first term in equation (17) 
is the expected noise power ((Pat) = 4 because we have summed four power spectra each 
normalized to unity), while the second term is the expected signal power. Here, we have 
assumed that the signal waveform is sinusoidal. Ng and Ni are determined for each data 
set before beginning the pulsation search. Therefore, the sensitivity of each search is 
determined as a limit on the pulsed fraction T . 

We now calculate the detection sensitivities for the eight data sets searched. For a 
given observation, the detection sensitivity is expressed as the minimum pulsed fraction 
required to produce a spectral power exceeding Pdet with high confidence. Instead of simply 
solving equation (17) for T using the detection threshold power P^et = 43, we must allow 
for statistical variation of the spectral power produced by a source of a given strength, as 
well as variation in the recovered power due to the discrete nature of the search trials. 

To determine the detection sensitivity, we must consider the probability distribution 
of power in a spectral bin containing a signal plus noise. If, in the absence of noise, the 
signal produces a power Pgig, then the power P in the presence of noise will be distributed 
according to 

oo TO+n— 1 

p„(P>Po;Psig)=exp[-(Po + P,ig)]5^ P^P^J{k\w>) (18) 

m=0 fc=0 

(Groth 1975), where n is the number of independent bins that were summed to produce P 
(i.e., in our case n = S = 4). Thus, Pn{P > Po! Psig) is the probability that the power will 
exceed Pq in the spectral bin containing the signal. It is important to note that noise power 
and signal power are not simply additive; P will not always exceed Pgig (see Vaughan et al. 
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(1994) for a discussion of this point). For a given Ng and Nt, we would like to find the 
pulsed fraction T that is 95% likely to produce power P > P^et- 

To ensure that our reported sensitivity is as conservative as possible, we will not 
naively invert equation 18 for Pgig with Pq = Pdct- Instead we consider the worst case 
scenario covered by our search parameter space. Due to the discrete binning of photons 
in the construction of the time series, the recovered signal power in the FFT falls off with 
increasing frequency. A signal near our maximum search frequency — 1024 Hz, which is 
half the Nyquist frequency of our power spectra — will produce 81% of the power that 
would be recovered from a low frequency signal of the same intrinsic strength. Also, due to 
our discrete grid of acceleration trials and frequency trials, we are only 98% likely to recover 
more than 77% of a signal's available power. Thus, in our worst case, we are 98% hkely to 
recover at least (0.77)(0.81)Psig = 0.62Psig. 

We account for this reduction in recovered signal power approximately by solving 
equation 18 in the form 

p„(P>Pdet;0.62P,ig)=0.97. (19) 

The probability of our recovered signal power exceeding 0.62Psig and our total spectral 
power P exceeding P^et is then approximately (0.98) (0.97) ft! 0.95, giving our desired 95% 
confidence. 

Solving equation 19 (numerically) results in Pgig = 93.6. For a given observation, the 
minimum required pulsed fraction is then given by ^ = {ANtPsig/NgY^"^ and converted to 
RMS by dividing by 2^/2. Thus, our 95% confidence detection sensitivities for observations 
8a, 10a, 10b, 11a, lib, 11c, 12a, and 12b are given by 0.60%, 1.2%, 1.2%, 2.7%, 2.9%, 3.0%, 
5.6%, and 10% respectively. We verified these numbers with Monte Carlo simulations. 

Our search included spin frequencies up to 1024 Hz. Had we restricted our search 
to a smaller range of frequencies, corresponding to the previously reported pulsations. 
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our detection sensitivities would not have been drastically different. To achieve the same 
detection significance, with far fewer trials, we require only a slightly reduced spectral 
power. And since the pulsed fraction limit depends essentially on the square root of 
the required power, the detection sensitivities are not terribly sensitive to the number 
of frequency trials in the search. Given Fox's best estimate for the pulsation frequency 
549.76^q'q3 Hz (D. Fox 1999, private communication), if wc had confined our search to a Scr 
range (about 549.76 Hz and 274.88 Hz), our detection sensitivity in observation 10a would 
have improved to 1.0%, while in observation 12b the hmit would have been 8.8%. 

2.5. Verification of Procedures 

Simulated data were used to test our search codes and to verify the sensitivity of the 
search. As an additional check, we applied our acceleration method to SAX J1808. 4-3658 
— the accreting millisecond pulsar (Wijnands & Van der Klis 1998; Chakrabarty & Morgan 
1998). Since our search was for a similar pulsar, with ~ few hundred Hz frequency in a 
~ few hour orbit about a low-mass companion, SAX J1808. 4-3658 provides a useful test of 
our acceleration method. 

We analyzed an RXTE/PCA observation (observation #30411-01-01-04). We selected 
energy PHA bins corresponding to 0.4-17 keV from data of the same type we use in our 
search (E125). We barycenter corrected the arrival times with the position given by Roche 
et al. (1998), and used 2048 seconds of data beginning at Mission Elapsed Time (MET, 
seconds since 1994.0) 135395082. 

The pulsations in these data are easily detected — the ~ 401 Hz signal is obvious in an 
unacceleratcd PDS, with significant power detected in each of about 30 adjacent frequency 
bins, using 2048 seconds of data (Figure 1, dotted line). The highest single bin power 
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represents about 22% of the total signal power in the extended feature. 

Using the known orbital parameters (Chakrabarty & Morgan 1998) to remove 
the acceleration from the data, the spreading of the signal power is reduced to about 6 
independent Fourier bins around the correct rest frequency of the pulsar, 400.975 210 6(8) Hz, 
with ~ 85% of the total power in a single spectral bin (Figure 1, solid line). Oversampling 
reveals that the spreading that remains is consistent with the expected sinc-function 
response of the (discrete) FFT to a signal with the pulsar's rest frequency (Van der Khs 
1989). Thus, we confirm that our technique successfully removes the effect of the doppler 
shift due to the pulsar's orbital motion. 



2.6. Seeirch Results 

No candidate pulsations from any of the Aql X-1 observations exceeded our 
predetermined detection threshold. Using our largest detected search power, we can 
calculate upper limits on the strength of any pulsar signal contained in the various data 
sets searched (Vaughan et al. 1994), provided its orbit and spin frequency are covered by 
our search phase space. These upper limits are more restrictive than the a priori detection 
sensitivities. 

The upper limit calculation is essentially the same as the detection sensitivity 
calculation, with the detection threshold power P^ct replaced by the maximum observed 
power Pmax = 30.96. The results are shown in Table 3. For example, with 95% confidence, 
we can say that there is no sinusoidal pulsar signal in observation 10a with a fractional 
RMS amplitude greater than 1.0%. The limits we can place on the fractional RMS 
become less stringent as the source fiux decreases, increasing to an upper limit of 9.0% 
during observation 12b, when Aql X-1 was observed to be faintest (~1 mCrab; 13 000 
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PCA counts s-^=l Crab). 

3. DISCUSSION 

We find no evidence of pulsations from Aql X-1 as it fades into quiescence. The upper 
limits on the pulsed fraction for the eight data sets searched range from 0.52% to 9.0% 
RMS. If, indeed, during the period identified as the "propeller phase" (Campana et al. 
1998; Zhang et al. 1998a) the NS magnetic field significantly modifies the accretion flow- 
geometry in the vicinity of the NS, the apparent absence of pulsations does not support the 
hypothesis that the quiescent emission is due to continued accretion. Our results cannot 
completely rule out such accretion, however, since it is possible that the geometry of the 
system may not lead to detectable pulsations; for example, the magnetic axis may be very 
nearly aligned with the rotation axis, or the rotation axis may point directly towards the 
earth. The non-detection of pulsations coupled with the observation of a thermal spectrum 
during quiescence (Rutledge et al. 1999) favors the interpretation that the quiescent 
luminosity is not due to accretion, but rather to a hot NS core (Brown et al. 1998). 

We thank Stuart B. Anderson and Thomas A. Prince for helpful discussions. We also 
thank Lars Bildsten, Sergio Campana, and an anonymous referee for valuable comments on 
the manuscript. This work was supported by NASA Grant #NAG5-3239. 
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Fig. 1. — SAX J1808. 4-3658. The dotted line shows a section of the unaccelerated power 
spectrum. The spreading of the signal power due to the orbital doppler shift is evident. 
The solid line shows the power spectrum after removal of the doppler shift using our ac- 
celeration technique. The signal power has not been confined to a single bin because the 
pulsar frequency lies between two discrete Fourier frequencies. During this observation, 
SAX J1808. 4-3658 passes through longitude = 0. The peak in the uncorrected spectrum is 
therefore maximally offset from the pulsar's rest frequency, A/max = Pfo = 0.022 Hz. 
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Table 1: RXTE Aql X-1 Analyzed Observations 



Observation 


Start Time / End Time 


Start Time 


Start Time 




(UT) 


(JD) 


(MET) 


20098-03-08-00 


1997-03-01 21:33 1997-03-02 00:21 


2450509.398 


99869572 


20098-03-10-00 


1997-03-05 22:55 1997-03-06 01:28 


2450513.456 


100220227 


20098-03-11-00 


1997-03-08 21:31 1997-03-09 01:11 


2450516.394 


100474022 


20098-03-12-00 


1997-03-10 21:25 1997-03-10 23:55 


2450518.389 


100646431 
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Table 2: Searched Parameter Space 



Parameter Parameter Space 

Porb 18.91 - 19.45 hr 

vj^s sin i — 130 km 

Initial Orbital Phase 0—1 cycles 

Pulsar Spin Frequency 0.5 — 1024 Hz 
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Table 3: Observation Parameters and Search Results 





iVill/ ± 


ooui ce 


ijdCKgi ounci 


ri actional rtivio 


Observation 


Analyzed 


*^ount riate 


v^ount riate 


Ampiitucie 


Sc. 


yyooyo / z-yso / idzu 


«7n R-un s 

D/U.D±U.o 


Q7 Q 


< u.oz/o 




iUUZzUzz (-iUUZZZZ ( 


iOO.ZztU.D 


Qc: 
OO.o 


< i.U/o 


iUD 


lUUZZDZZ / -lUUZZoZ / 


1 77 /( 4_n p. 


OO.o 


< i . i /o 


11a 


100474022-100476070 


52.6±0.5 


35.7 


< 2.3% 


lib 


100479522-100481570 


48.6±0.5 


35.7 


< 2.5% 


11c 


100485082-100487130 


45.5±0.5 


35.7 


< 2.6% 


12a 


100646431-100648479 


19.0±0.5 


31.2 


< 4.9% 


12b 


100652431-100654479 


9.2±0.5 


31.2 


< 9.0% 



"counts (13 000 counts s'^ = 1 Crab) 
''Systematic uncertainty ~ 0.5 counts s~-^ 
'^Upper limits are 95% confidence 



